Robustness of Quadratic Solitons with Periodic Gain 
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We address the robustness of quadratic solitons with periodic non-conservative perturbations. 
We find the evolution equations for guiding-center solitons under conditions for second-harmonic 
generation in the presence of periodic multi-band loss and gain. Under proper conditions, a robust 
guiding-center soliton formation is revealed. 



Multicolor optical soliton formation mediated by cas- 
cading of quadratic nonlinearities has been demonstrated 
experimentally during the last few years in a variety of ge- 
ometries and frequency-mixing processes, in settings for 
spatial, temporal and spatio-temporal trapping of light. ^ 
Soliton signals exist in particular in the process of second- 
harmonic generation (SHG) that is addressed here, where 
solitons form in waveguides and in bulk crystals by the 
mutual trapping between the fundamental frequency and 
second-harmonic waves. Multidimensional soliton fami- 
lies exist above a threshold light intensity for all values of 
the phase-mismatch between the waves, and most of such 
solitons have been shown to be dynamically stable under 
propagation with the equations that model the ideal light 
evolution under conditions of focused and pulsed SHG. 
Adiabatic soliton decay and amplification in the presence 
of weak loss or gain, have been also studied. ^""^ 

In this Letter we address the robustness of quadratic 
solitons against strong, but periodic non-conservative 
perturbations. To start the program we consider soli- 
ton formation in the presence of multi-frequency losses 
and large, but rapidly-varying periodic gain. Our goal 
is to derive the corresponding guiding-center evolution 
equations and to expose the robustness of the existing 
solitons under proper conditions. We believe that the re- 
sults reported bear a generic fundamental interest to the 
robustness of quadratic soliton formation in structures 
with periodic non-conservative perturbations. Moreover, 
they might find direct applications in reduced models of 
multi-color laser systems with intracavity frequency gen- 
eration, including self- frequency doubling schemes, oper- 
ating in the solitonic regime. 

Here the focus is on solitons formed in one-dimensional 
structures under conditions for non-critical type I SHG, 
but the analysis can be extended to different physical set- 
tings. The evolution of the slowly-varying envelopes of 
the light waves in the presence of multi-frequency band 
loss and periodic gain can be described by the reduced 
equations'' 



.dai ai d^ai 
.da2 a2 9^a2 



■ al 02 exp(-i/5^) = iri(^)ai, (1) 
a? exp(i/?0 =*r2(0a2, (2) 



where ai and 02 are the normalized amplitudes of the 
fundamental frequency (FF) and second-harmonic (SH) 
waves. In the case of spatial solitons, ai = —1, and 
Q!2 = —ki/k2 — —0.5, where k^, with v — 1,2 are the 
linear wave numbers at both frequencies. In the case 
of temporal solitons, stand for the ratio between the 
group- velocity-dispersions existing at both frequencies. 
The transverse and longitudinal coordinates are normal- 
ized to the beam or pulse width and to the diffraction 
or dispersion length at the fundamental frequency, re- 
spectively. The parameter (3 is the scaled phase mis- 
match, and ^uiO stand for periodic gain and loss. Let 

riy(C) — ri''^ + Ri^iOi where r° are the average gain or 
loss at the FF and SH frequencies, and R,y{£,) are periodic 
functions, with period L <^ 1, and zero mean. 

To derive the evolution equations of the guiding-center 
solitons in the presence of the periodic gain, we first use 
the approach originally introduced by MoUenauer, Evan- 
gelides and Haus for the case of Kerr solitons propagating 
in optical fibers. * Let Ai = ai and A2 = 02 exp(— 
The explicit periodic gain can be removed from equa- 
tions (l)-(2) by making the transformations A^{^,s) — 
U^{^,s)e^p[J^ R^{C)dC]. Substitution in (l)-(2) leads 
to the evolution equations 
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(3) 



zrf C/2-/3C/2 + e3(e)f/f = 0, (4) 



where the resulting longitudinally-varying nonlinear co- 
efficients write 



£2(0 = exp 
£3(0 = exp 



i?2(^')dr 
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(5) 
(6) 



To proceed further one now assumes that the wave evo- 
lution over a period of the map is mostly dictated by the 
gain and loss, which are responsible for fast amplitude 
oscillations of the fields A^, in addition to the residual 
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effects induced by the nonlinearity. Therefore, Ui, vary 
slowly over a period of the map, so that one can average 
equations (3)- (4) to approximately get 



expand ( jisj ) in powers of the corresponding parameters 
to get, at order 



.dUi aid^Ui (0) 

.dU2 a2d^U2 .„(0),, r,rr I tt2 n 
*^"T^" 2 ^2-/3C/2 + 773 C/l -0, 



where 



772,3 =< £2,3(0 >= 



e2,3(e)dO 



(7) 
(8) 

(9) 



are the averaged nonlinear coefficients over a period of 
the map. Therefore, the evolution equations for the 
slowly- varying averaged fields U,y,Q{^,s), defined as 



C/,,o(Os)=<exp[/ i?,(0)dO] >C/.(e,s), (10) 
Jo 



write 
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(11) 



in"U2fi- PU2fi + vUto-0, (12) 



dS, 2 c)s2 
where 77 = ri2r/3/rjl, with ryi =< ei(0 >i and 



ei(0 = exp 



(13) 



One thus concludes that under the conditions where 
the approximations used to derive (ll|)-(p^) hold, the 
guiding-center solitons are given by those existing with- 
out gain and loss when F^P^ = 0,^ or otherwise with com- 
pensated gain and loss,^*^ but with properly renormal- 
ized amplitudes. To expose the value and dependencies 
of the renormalization coefRcient rj, consider the illustra- 
tive map i?,/(0 — Gy H{£^), where H{^) is the Heaviside 
function 



H{0 



-1, nL < ^ < {ji + l/2)L, 



1, 



(n + l/2)L<C< {n + l)L. 



(14) 



Substitution in the above expressions and performing all 
the averages yields 



where 



F[G2L/2]F[{2Gi-G2)L/2] 
F2[GiL/2] 

1 - cxp(-/^) 



m = 



(15) 



(16) 



Of particular interest are the weak maps with /i ^ 1, 
corresponding to G^L <C 1. In such cases, one might 



(Gi - G2 )^L^ 
48 



ry ~ 1 + 



(17) 



Note that when Gi = G2, equation (|15|) gives r/ = 1, 
regardless the value of the map period L. However, this 
does not mean that the guiding-center evolution is given 
by (|rTl)-(|l^) with 77 = 1 at all orders of L, because 
the derivation of (pT|)-([l^), hence the governing equa- 
tions themselves, arc only intended to hold when i <C 1. 
To elucidate the applicability limits of (pl])-(p^, it is 
worth deriving the guiding-center evolution equations 
using more mathematically systematic approaches. ^^'^^ 
Next we outline the outcome of the asymptotic expansion 
method developed by Kivshar and co-workers for simi- 
lar problems but for Kerr solitons. Such approach has 
been recently employed to obtain guiding-center evolu- 
tion equations of light signals propagating in quasi-phase- 
matched quadratic structures. ^'^ 

Assuming perfect periodicity of the functions Ri, {^) as 
above, one can express the periodic gain and all fields as 
Fourier series of the form 



■Ri^(C) = y^d,y,Kexp(wAg), 



^1^(0 s) = X! ^''■"(0 s) exp(mAO, 



(18) 



(19) 



where A = 2tt/L and the Foruier coefficients U^,n{^,s) 
vary much slower than the corresponding carrier 
exp(mAO- For the map (|l^ ) considered here one has 



diy.2n — 0, 



Ji/,(2n+l) 



2G, 



i7r(277+ 1)' 



(20) 



Assuming the amplification maps to vary rapidly over a 
diffraction/dispersion lenght, i.e. L <^1, so that the spa- 
tial frequency A is large enough, one can expand the har- 
monic amplitudes C/i/,„^o, as power series of 1/A. Namely, 



(m) 



m>0 



(21) 



Substitution in the governing equations and matching 
leading-order contributions leads to 



^(1) r| 



(22) 



The guiding-center evolution is obtained by solving re- 
cursively for the higher-order contributions i^i.n''^'' and 
substituting the corresponding expressions in the evolu- 
tion equations for the average fields C/^.o- At order l/A^, 
one exactly finds (|ll-(|l2|), with (0). 

Equations ( pT| ) , ([12|) and (|l^) are the central result of 
this paper. First, they reveal that when Gi = G2, at 
order the guiding-center evolution equations are iden- 
tical to those without the periodic gain, regardless the 
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value of the amplitudes Gj^. This is a remarkable result, 
that emphasizes the robustness of quadratic solitons with 
the type of perturbations considered. Second, the typical 
values of the gain coefficient that is obtained in Erbium- 
doped lithium niobate or potassium titanyl phosphate 
around the third telecommunication window centered at 
1.55 /im fall in the range 0-2 dB/cm.^'-'^^ In the case of 
spatial solitons, this yields values of the normalized gain 
coefficient of the order of G^, ~ 10~^.^ With such val- 
ues and letting L ~ 10~^, one always obtains negligible 
corrections of the order of [77 — 1] ~ 10~^. 




FIG. 1. Typical evolution of guiding-center quadratic soli- 
tons in structures with rapidly-varying periodic gain. In 
(a)-(b): The input is a renormalized stationary soliton so- 
lution. In (c)-(d): The input contains only fundamental fre- 
quency light with a Gaussian shape, (a), (c) : Evolution of 
the FF; (b), (d): Evolution of the SH. Conditions: /3 = 3, 
/ = 30, Gi = 20, G2 = 0, L = 0.2. 

To confirm that under proper conditions (pT|)-(p^) 
hold, we solved (l)-(2) numerically for different input 
conditions and maps. Figure 1 shows typical exam- 
ples of the outcome, when r[,°^ = 0. The plots corre- 
spond to the phase-mismatch /3 = 3, but analogous re- 
sults were obtained for other values. To emphasize that 
guiding-center solitons form with gigantic values of the 
gain-loss amplitude, provided that the period of the map 
is small enough so that the guiding-center approach is 
justified, we display results for a map with Gi = 20, 
G2 = and L = 0.2. Figures l(a)-(b) show the prop- 
agation of a guiding-center soliton excited by the input 
A^(0) = A^{0)/p^, where Aj^(O) are the corresponding 
stationary solitons existing without gain and loss with 
energy flow I = ^ + = /(liip + |i2|')ds = 30,^ 
while pi = (772773)^^^, and £2 = 772, are the renormal- 
ization factors dictated by For the map consid- 

ered one has pi ~ 0.5, so that in the renormalized input 
the FF energy is enhanced by a factor of four. Figures 
1(c)- (d) show the excitation of a guiding-center soliton 
in the same map but with only FF input light carrying 



the same energy flow as above in the Gaussian shape 
Alio) = (2/^)1/4 exp(-s2). Here h = h/pl 
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Propagatbn Coordinate 

FIG. 2. Detail of the evolution of the peak amplitudes of 
the signals displayed in Figs. l(a)-(b), over the first five peri- 
ods of the map. 

Figure 2 shows the detailed evolution of the field am- 
plitudes in the case shown in Figs. l(a)-(b). Because of 
the large excursions of the amplitudes, the actual evo- 
lution differs slightly from that predicted by equations 
(pl])-(p^). Such departures are responsible, e.g., of the 
small resonance peaks appearing in the Fourier spectra 
of the FF and SH evolving signals, as shown in Figure 3. 
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FIG. 3. Fourier spectra of the FF and SH signals displayed 
in Figs. l(a)-(b), at ^ = 40. 

When r^°^ 7^ 0, but the net loss at one frequency is 
compensated for by the presence of a net gain at the 
other frequency band, the system (pT|)-(p^ allows sta- 
tionary soliton solutions,^" which are chirped. However, 
an interesting result revealed by the guiding-center evo- 
lution equations ([Tl|)-([l^) is that, in the absence of a 
net gain or loss (i.e., T)^' = 0), on average the guiding- 
center quadratic solitons are chirpless. Accordingly, in 
contrast to what is found with solitons in other periodic 
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systems, e.g., with dispersion-managed solitons propa- 
gating in optical fibers, one concludes that with chirp- 
less input conditions the guiding-center quadratic soli- 
tons with periodic gain are best excited when the first 
domain has the whole nominal length. The outcome of 
our numerical simulations confirms that such is indeed 
the case, as shown in Figure 4. The plot shows the evo- 
lution over a period of the map of the quantities 



cm 
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/s2|A^|2ds 



(23) 



where ^ stands for dAy/ds. The evolution displayed 
corresponds to a guiding-center soliton similar to that 
shown in Figs. l(a)-(b), but for Gi = G2 = 10. 
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FIG. 4. Typical evolution of the integral chirp of the guid- 
ing-center solitons over a period of the map. The plot shows 
the evolution from ^ = 19.8 to ^ = 20. Conditions: /3 = 3, 
/ = 30, Gi = 10, G2 = 10, L = 0.2. The input is the corre- 
sponding renormalized stationary soliton solution. 




FIG. 5. Unstable evolution when the guiding-center ap- 
proach fails. The input is a renormalized stationary soliton 
solution. Same conditions as in Figs. l(a)-(b), except that 
here Gi = 2 and L = 1. (a): Evolution of the FF; (b): Evo- 
lution of the SH. 

Analogous results than those shown above were ob- 
tained for a variety of values of the phase-mismatch /3, 
the input energy flow /, the signal shape, and the map 
amplitude. Naturally, this is so provided that the input 
conditions, the gain amplitude, and the period L make 



a guiding-center approach justified. ^^"^^ Otherwise, e.g., 
when L, 0^ or the input / are too large, the guiding- 
center evolution fails, yielding a totally new scenario. In 
particular, under such conditions resonance phenomena 
that make the wave propagation unstable can occur. Fig- 
ure 5 shows a typical example. The plot shows the un- 
stable evolution of the renormalized solution with / = 30 
in a map with Gi = 2, with L = 1. 

In conclusion, the evolution equations for guiding- 
center quadratic solitons propagating in structures with 
multi-frequency losses and rapidly-varying periodic gain 
have been presented. Under proper conditions, robust 
multicolor quadratic soliton formation has been revealed. 
Results might find applications to reduced models of mul- 
ticolor laser systems with intracavity frequency genera- 
tion, including self- frequency doubling schemes, operat- 
ing in the solitonic regime. Extension of the analysis to 
general maps that include the details of the laser struc- 
tures, including the pump light-matter interaction,^'-'^'* is 
worth investigating. 
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